Observation of B+ ii:i(1270)+7 
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We report the observation of the radiative decay — + Ki (1270)^7 using a data sample of 
140 fb"^ taken at the T{4S) resonance with the Belle detector at the KEKB e'^e collider. We 
find the branching fraction to be B{B+ -> iCi(1270)+7) = (4.3 ± 0.9(stat.) ± 0.9(syst.)) x 10"^ 
with a significance of 7.3a. We find no significant signal for B^ — > J('i(1400)^7 and set an upper 
limit B{B+ A'i(1400) + 7) < 1.5 x lO"'^ at the 90% confidence level. We also measure inclusive 
branching fractions for B^ K^Tv^Tv^'y and B" — > K^TT^n^'y in the mass range 1 GeV/c^ < 
M^+(o,,+ ^- < 2 GeV/c^ 
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Radiative B decays that occur through the flavor 
changing neutral current process b ^ sj have been one 
of the most sensitive probes in the search for physics 
beyond the Standard Model (SM). The first observed ex- 
clusive radiative decay mode was B K*j Q, 0- The 
second mode was B i4r|(1430)7, for which evidence 
was reported by CLEO and Belle [jj]. Theoretical pre- 
dictions for the branching fractions of the unobserved 
exclusive decays can be found in Refs. 0,0. The modes 
B i^i (1270)7 and B i^i (1400)7 (^i ^ ^tttt) 
can be used to measure the photon helicity, which may 
differ from the SM prediction in some models beyond 
the SM 0. The neutral mode B° i5:i(1270)°7, 
A'i(1270)° ^sP^ would also be useful to measure 
time-dependent CP violation that may arise from new 
physics 0. 

In this paper, we report the observation of B^ 
ifi(1270)+7, which is the first radiative B meson de- 
cay mode that involves an axial-vector resonance. We 
study radiative decays in the K~^tt^tt^j and Kgir^ir^j 
final states, where we search for resonant structure in 
the K'!r'^Tr~ system 0- We also report inclusive mea- 
surements of B^ i^+7r+7r~7 and K^7T~^Tr~j, 
and the results of a search for i?+ i<ri(1400)+7. The 
analysis is based on a data sample of 140 fb~^ (152 mil- 
lion BB pairs) taken at the T(4S') resonance with the 
Belle detector at the KEKB e+e-colHder 

The Belle detector consists of a three-layer silicon ver- 
tex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cerenkov counters 
(ACC), a barrel-like arrangement of time-of-flight scintil- 



lation counters (TOF), and an electromagnetic calorime- 
ter comprised of CsI(Tl) crystals (ECL) located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An instrumented iron flux-return (KLM) for 
K1 and muon identification is located outside of the coil. 
The detector is described in detail elsewhere 0. 

The photon candidate is the highest energy photon 
cluster measured with the barrel ECL (33° < 6^ < 128° 
in the laboratory frame). In order to reduce the back- 
ground from TT^/ry — > 77 decays, we combine the pho- 
ton candidate with all other photon clusters in the event 
with energy greater than 30 MeV (200 MeV) and reject 
the event if the invariant mass of any pair is within 
±18 MeV/c^ (±32 MeV/c^) of the nominal 7r° (77) mass. 
These correspond to ±3(7 windows, where a is the mass 
resolution. We refer to this requirement as the Tr'^/ry veto. 

Charged tracks are required to have momentum in 
the center-of-mass (cm.) frame greater than 200 MeV/c 
and to have an impact parameter relative to the inter- 
action point of less than 5 cm along the positron beam 
axis and less than 0.5 cm in the transverse plane. The 
charged kaon candidate is identified using a likelihood ra- 
tio combining the information from the ACC, TOF, and 
CDC sub-detectors; the remaining charged particles in 
the event are used as pion candidates, unless the track 
has been identified as an electron, muon, or proton. 

For neutral kaons, we use tt+tt^ candidates that 

have invariant masses within ±30 MeV/c^ of the nominal 
K'g mass and a cm. momentum greater than 200 MeV/c. 
The two pions are required to have a common vertex that 
is displaced from the interaction point. The momen- 
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turn direction is also required to be consistent with the 
Kg flight direction. 

In order to study B Kij {Ki s- Kmr), we 
first reconstruct B —>■ Ktt^tt^j inclusively, without 
any requirement for the structure of the Ktt'^tt' sys- 
tem. We select A'tt+tt^ combinations in the mass 
range 1 GeV/c^ < Mk^+^- < 2 GeV/c^. Given the 
i^TTT+TT" system and the photon candidate, we identify 
B meson candidates using two independent kinematic 
variables: the beam-energy constrained mass Mbc = 



lc?Y ~ \'PbI'^V ^''^d the energy difference Ai? = 
^^boam' where El^^^ is the beam energy 
and 'Pq is the momentum of the B candidate in the 
cm. frame [loj . The B momentum is calculated as 

Vb = V^^. + 1^ X (^bcam - ^K.-n) ^ Order to im- 
prove the A/bc resolution. We select B candidates within 
-0.1 GeV < A£; < 0.08 GeV and Mbc > 5.2 GeV/c^. If 
there exist multiple candidates, we choose the candidate 
with the highest confidence level for the if "•"tt+tt" vertex 
fit (7r"'"7r~ vertex in the neutral case). 

The dominant background comes from hadronic con- 
tinuum (e+e~-^ gg, q — u,d,s,c). To suppress this 
background, we use two variables: the B flight direc- 
tion (cos^Ij) and a Fisher discriminant ^3 built from a 
set of shape variables For signal, the B flight direc- 
tion follows a 1 — cos^ 9'^ distribution while that of qq is 
nearly uniform. The likelihood function jOf^g^^^ is mod- 
eled as a 2"^^ (1*^*) order polynomial for the signal (con- 
tinuum background) from MC samples. For the shape 
variables, we use 16 modified Fox- Wolfram moments [l3| 
calculated for the following groups of particles: 1) par- 
ticles that form the signal candidate, 2) the remaining 
charged particles, 3) the remaining neutral particles, and 
4) a hypothetical particle for the missing momentum of 
the event. The Fisher discriminant is obtained from these 
moments and the scalar sum of the transverse momen- 
tum. The likelihood function -Cl'/^"' is modeled as a 
bifurcated Gaussian function both for the signal and the 
continuum background from MC samples. 

These likelihood functions are then combined to form 

77 /'COS /'Fisher / / /"Cos ^pjgj^cj. /'Cos 0^ ^pjgbcr\ 

/Cs — i-s '-S I \'-S '-B '"B )■ 

We determine the TZs requirement by maximizing 
I \J^s + ^B , where Ns and Nb are the expected 
number of the signal and background events, respec- 
tively, in Mbc > 5.27 GeV/c^. For this purpose, we 
use B ^ ii:i(1270)7 and B i4:i(1400)7 signal Monte 
Carlo (MC) simulated data, assuming all the B A'17 
branching fractions are 1 x 10~^. We find TZs > 0.9 is 
the optimal requirement. This requirement retains 47% 
of the signal events while rejecting 98% of the continuum 
background events. 

The signal yields for B K-k^tt^^ are extracted from 
a binned maximum likelihood fit to the Afbc distribu- 
tion. In addition to the continuum background, we con- 
sider four B decay background sources: known B decays 



through the b ^ c transition (referred to as the b c 
background), hadronic B decays through the b u, 
d or s transitions (charmless background), B K*j 
background, and radiative 6 — *■ s decays to final states 
other than K*j and K'!T^Tr~^ (other b back- 
ground). To suppress the B —t K*"f background, we 
reject the event if AE and A/bc calculated from either 
Kn-/ combination satisfy -0.2 GeV < A£^ < 0.1 GeV 
and A/bc > 5.27 GeV/c^. The signal A/bc distributions 
are each modeled as a Gaussian function; its width is 
fixed using a data sample oi B D{-^ K'k-k)'k decays, 
treating the primary pion as a high energy photon. The 
shapes of the background A/bc distributions are deter- 
mined using large MC samples. We find that the sum of 
the continuum and b ^ c backgrounds is described by an 
ARGUS function , a smooth functional form that has 
a kinematic threshold at half of the center of mass energy. 
Charmless decays, B A'*7 and other b sj decays 
are modeled as a sum of an ARGUS function and a Gaus- 
sian function. The normalization of the continuum plus 
b ^ c background is floated in the fit; the normalization 
of the other three components are fixed in the fit. 

The fit result is shown in Fig. ^ For the /3+ 
/ir+7r+7r^7 mode, we obtain 318 ± 22 events with a sig- 
nificance of 16(T, where the significance is defined as 
\/—2 ln(£o / 'Cmax) , and £max and £0 denote the maxi- 
mum likelihoods of the fit with and without the signal 
component, respectively, and the significance includes 
systematic error. Similarly, we obtain 67±10 events with 
a significance of 8.3(t for the K^tt^tt^j mode. 

The systematic uncertainty related to the fitting pro- 
cedure is estimated in the following way. We vary the 
width and the mean of the signal Gaussian by the error 
of the B Dtt calibration sample. We vary the AR- 
GUS parameter of the continuum plus b ^ c background 
by the errors from fits to the MC sample and to a data 
sideband region defined as 0.1 GeV < AE < 0.5 GeV, 
then we take the quadratic sum of those errors. The 
B —>■ K*^ component is varied by the branching frac- 
tion uncertainty. The normalization of the other b ^ 
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FIG. 1: A/bc distributions for B+ K+Tr+-K~-f (left) and 
tt'^-k 7 (right). Curves show the continuum plus 
b c background component (dot-dashed), total background 
(dotted) and the total fit result (solid). 
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background component is varied within its respective un- 
certainty, estimated from the uncertainties in the to- 
tal 6 ^ S7 branching fraction [la and the fraction of 
if7r+7r~7 in the 57 final state [iQ]. For the charmless 
background we vary the normalization by ±100%. We 
also assign the uncertainty due to a possible fit bias as 
the error of the fit to the signal MC sample. The total 
fitting errors are 5.3% (12%) for the B+ K+tt+tt-j 
{B° K^TT+TT-j) mode. 

In order to decompose intermediate resonances that 
may be involved in the K~^Tr~^Tr~ final state, we per- 
form an unbinned maximum likelihood fit to the Mbc and 
^K+T^+n- distributions of the 5+ /ir+7r+7r~7 can- 
didates. There are many possible resonances that can 
contribute: A:i(1270), A:i(1400), A:2*(1430), A'*(1410), 
iir*(1680), and so on. We consider the first three reso- 
nances, and include an additional non-resonant B~^ 
K~^Tr~^Tr~j component. The B if^ (1430)7 compo- 
nent, which is already measured, is fixed in the fit. We 
model the Mjf+^+^- distribution of the ii'i(1270) res- 
onance as a sum of three decay chains, A'i(1270)+ 
K+p°, p° ^ TT+TT"; A:i(1270)+ ^ K*°Tr+, K*° 
K+TT-- and A:i(1270)+ -> Xo*(1430)°7r+, A'o*(1430)° 
K^n^ . The Mx+-K+-K- distribution for each decay chain 
is described by convolving the two relativistic Breit- 
Wigner functions of the resonances in the decay chain. 
The Mx+TT+Ti- distribution of the Ki{lAQ{)) resonance is 
modeled with a single decay chain, iiri(1400)+ /i*"7r+, 
AT*" AT+TT". The Mx+Tt+Tt- distributions of other 
6 S7, non-resonant B^ t:~ ^ and continuum 

plus b —>■ c backgrounds are modeled using the function 
(Po + Pix) exp{p2 +P3X + pix'^), where x = M^+^+^r-, 
and Pi {i — ... 4) are empirical parameters that are 
determined from MC samples. 

In order to enhance the i(ri(1270) component, we se- 
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FIG. 2: Mbc distribution for 1.2 GeV/c^ < Mk+^+^- < 
1.4 GeV/c^ (left) and Mj^+^+^- distribution for Mbc > 
5.27 GeV/c^ (right) of the /«'i(1270)+7 enriched sample with 
0.6 GeV/c^ < < 0.9 GeV/c^ Curves show the projec- 

tions of the fit results for the continuum plus b ^ c back- 
ground component (dot-dashed), total background without 
and with the /5r|(1430)'''7 component (dotted and dashed, 
respectively), A'i(1270)+7 (thin solid line) and Ai(1400) + 7 
(hatched) components, and the sum of all components (thick 
solid line). 



lect events with tt+tt mass in the p mass region, 
0.6 GeV/c^ < M^^ < 0.9 GeV/c^ (i3(A:i(1270) -> Kp) 
= (42 ± 6)% being much larger than B(A:i(1400) -> Kp) 
= (3 ± 3)% 113). The fit result is shown in Fig. H 
We find 102 ± 22 events for the B+ A:i(1270)+7 
component with a significance of 7.3cr, while we fix the 
yields of the K*j, other b sj and charmless back- 
ground components to be 2.3, 10.2 and 2.0 events, re- 
spectively. Similarly, we select the events with the AT+tt" 
mass in the AT*" mass region, 0.8 GeV/c^ < Mktt < 
1.0 GeV/c^, in order to enhance the A'i(1400) compo- 
nent (6(A'i(1400) A'*7r) (94 ± 6)% being much 
larger than B{Ki{1270) K*n) = (16 ± 5)% O). We 
find 23 ± 26 events for the B+ -> A:i(1400)+7 compo- 
nent. Since the Ari(1400)+7 component is not signifi- 
cant, we set a 90% confidence level upper limit on the 
signal yield, iVgo, which is calculated from the relation 
J^''"' C{n)dn — 0.9 C{n)dn, where C{n) is the likeli- 
hood function with the signal yield fixed at n. 

The systematic uncertainty related to the fitting proce- 
dure includes the same contributions as for the K~^Tr~^Tr~ 
fit. The sub-branching fractions for the A'i(1270) and 
the branching fraction for the A'j (1430)7 ^-s well the un- 
certainties in these quantities are taken from Ref. 0| 
and included in the fitting procedure. We also assign 
the uncertainty due to incorrectly reconstructed signal 
events, which is determined from a MC study. The to- 
tal unbinned fitting errors are 20% (32%) for the 5+ 
A'i(1270)+7 {B+ A:i(1400)+7) mode. 

The reconstruction efficiency is obtained from the 
signal MC samples. For the inclusive measurement, 
we calculate the efficiency using the Ari(1270)"''7 and 
Ar+7r+7r~7 non-resonant MC samples, and weighting ac- 
cording to the measured ratio of yields in those modes 
(A:i(1400)+7 and A:^(1430)+7 contributions are ne- 
glected). We assume the same ratio for the B'^ 
A'°7r+7r^7 efficiency. We estimate the systematic er- 
rors due to photon detection (2.8%), tracking (1% per 
track), charged particle identification (0.5% per parti- 
cle). Kg reconstruction (4.5%) and the likelihood ratio 
and 7T°/r] veto (6.1% for A:+7r+7r-7, 3.2% for K^tt+tt-j). 
The total efficiency error is 7.6% (6.6%) for the B^ 
A'+7r+7r~7 {B^ A'°7r+7r~7) mode. 

Using the signal yield and efficiency we find 

BiB+ -> A:i(1270)+7) = (4.3 ± 0.9 ± 0.9) x 10"^ (1) 

where the first (second) error is statistical (systematic) 
assuming that the production of 5+ and 5° in T(4S') 
decays is equal. Wc also measure the inclusive branching 
fractions of S ^ A'7r+7r~7 given in Table HJ Wc find 
that B{B'^ — ^ Ar+7r+7r~7) is consistent with the previous 
measurement with a significantly improved error |l8l | ; the 
neutral one is measured with a similar branching fraction. 

Similarly, we perform an unbinned maximum likeli- 
hood fit to the B^ A'°7r+7r~7 candidates to decompose 
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TABLE I: Yields from fits, detection efficiencies, branching fractions with statistical and systematic errors, and significances. 
The 7^1(1270) and Ki{UOO) are reconstructed as A'tt+tt^ states which come through their sub-decay channels. 







yield 


efficiency 


branching fraction(B) 


significance 


B+ 


^ J<:i(1270) + 7 


102 ± 22 


(1.56 ±0.12)% 


(4.3 ±0.9 ±0.9) X lO-'' 


7.3cr 


B° 


^ J^i(1270)°7 


15.0^1',' 


(0.34 ± 0.02)% 


< 5.8 X 10"^ 




B+ 


^ A"i(1400) + 7 


23 ± 26 


(2.68 ± 0.20)% 


< 1.5 X 10"^ 






^ 7^i(1400)°7 


7 0+16.3 

-'•8-13.8 


(0.61 ± 0.04)% 


< 1.2 X 10"^ 




B+ 


—> K^n^n'^'y 


318 ±22 


(8.36 ± 0.64)% 


(2.50 ±0.18 ±0.22) X 10"^ 


16(7 






67± 10 


(1.82 ±0.12)% 


(2.4 ±0.4 ±0.3) X 10"^ 


8.3(T 



i4:i(1270)° and JCi(1400)°. Due to the limited statistics 
we do not obtain a significant result and set only upper 
limits. 

To conclude, we observe a new radiative decay mode, 
B+ A'i(1270)+7, with a branching fraction of (4.3 ± 
0.9(stat.) ± 0.9(syst.)) x 10~^, which is larger than the- 
ory predictions (0.5 ~ 2.0) x lO'^ 0,11. The rates for 
B Ki{l270)j and B /•i:i(1400)7 are sensitive to the 
magnitude and sign of the iiri(1270) — ifi(1400) mixing 
angle. The large rate for i?+ ifi(1270)+7 compared 
to B~^ — > ii'i(1400)+7 may be explained by a positive 
mixing angle This measurement of i? — *■ isTi (1270)7 
shows that in the future it will be possible to determine 
the photon helicity using B — *■ Ki^, Ki Kmr and that 
time-dependent CP violation using B^ iiri(1270)°7, 
ii'i(1270)° — *■ Kgp° decays can also be studied. We also 
measure similar branching fractions for the inclusive de- 
cays B+ if+7r+7r~7 and ^ K^tt+tt^j. The latter 
is measured for the first time. 
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